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ABSTRACT 

Signal reinforcements and cancellations due to standing waves may 
distort constant bandwidth combustion spectra. Cepstral techniques previ- 
ously applied to the ground-reflection echo problem are used to obtain 
smooth broad-band data and information on combustion noise propagation. 
Internal fluctuating pressure measurements made using a J47 combustor 
attached to a 6,44 m long duct are analyzed. Measurements made with Jet A 
and hydrogen fuels are compared. The acoustic power levels inferred from 
the measurements are presented for a range of low heat release rate operat- 
ing condHions near atmospheric pressure. For these cases, the variation 
with operating condition of the overall acoustic broad-band power level for 
both hydrogen and Jet A fuels is consistent with previous results showing it 
was proportional to the square of the heat release rate. However, the over- 
all acoustic broad-band power level generally is greater for hydrogen than 
for Jet A. 
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A duct area (0.0305m^) 

C(q) cepstrum 

Cq speed of sound, m/s 

Cp specific heat at constant pressure, 0/kg K 

f frequency, Hz 

G(f) rig system operation transfer function 

g(t) rig system operation function 

L duct length, m 

a location of microphone relative to duct exit, m 

(MW) gas molecular weight 

N(f) combustion process spectrum 

n(t) combustion process function 

P(f) pressure power spectrum 

p(t) pressure function. Pa 

p duct ambient pressure. Pa 

q quefrency, s 

gas constant, J/kg K 
T temperature, K 

t time, s 

W air mass flow rate, kg/s 

X r'uct axial distance, m 

Superscripts 

( * ] liftered function 
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INTRODUCTION 


As part of a combustior/ noise research program, pressure measurements 
were made in a combustion rig at the NASA Lewis Research Center using both 
Jet A and hydrogen fuels. At all test conditions pressure spectra showed 
the presence of organ pipe tones due to longitudinal duct resonance modes in 
addition to the broad-band combustion noise spectra. 

The presence of these organ pipe tones in the combustion spectra makes 
analyzing the experimental results difficult. For example, the spectra are 
not easily compared and the level of the integrated pressure may be deter- 
mined by the tones. Moreover, these tores are rig dependent. Consequently, 
in comparing experimental data from different rigs the tones must be 
extracted so that the general trend of the broad-band spectrum can be deter- 
mined. One approach to removing the tones is to har.d draw smooth lines 
through the experimental spectra (Ref. 1). These hand smoothed spectra can 
be effectively used to understand the parametric behavior of combustion 
noise; however, this is a tedious approach in an age of high speed computers. 

In this paper it is shown that cepstral techniques may be applied to 
obtain computer- smoothed combustion noise spectra. In addition, those 
techniques can be used to obtain information on sound propagation and 
reflection in the duct. This approach is possible because the measured 
signal can be considered to be composed of a combustion noise signal, 
reduced echoes of that signal due to reflections at the duct ends, and 
amplified echoes of that signal due to enhancement in the combustion zone. 
Each signal path is characterized by a different signal transfer function 
and a different signal arrival time. 
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The cepstral method is applied to both the hydrogen and Jet A combus- 
tion data. The major measured echo time delay is compared with the theo- 
retical echo time delay over a range of operating conditions for hydrogen 
and Jet A fuels. The computer-smoothed combustion duct pressure spectra are 
then presented for typical operating conditions. Last, the acoustic power 
levels calculated from the smooth spectra for the two fuels are compared. 


BACKGROUND 


The power cepstrum was first described by Bogart, Healy and Tukey as a 
technique to determine echo arrrival times (Ref. 2). The power cepstrum of 
a function is a Fourier transform of the logarithm of the power spectrum of 
that function. The cepstrum literature and application procedures are re- 
viewed by Childers, Skinner, and Kemerait in Ref. 3 and by Randall and Hee 
in Ref. 4. Cepstral techniques have been applied to the analysis of multi- 
cylinder engine exhaust noise (Ref. 5) and to the aeroacoustic ground 
reflection echo problem (Refs. 6 to 7). 

In order to accurately describe operations in cepstral analysis which 
might be confused with similar operations in spectrum analysis Bogart et al. 
(Ref, 2) introduced a now classical set of terms which are anagrams of the 
equivalent frequency domain terms. The ones used herein are: 

quef rency .frequency 

liftering filtering 

rahmonic. harmonic 

Note that the unit of quef rency is time. 
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In this report, cepstrum analysis is briefly discussed before the 
experimental investigation is presented. 


POWER CEPSTRUM ANALYSIS 


In the simple model used to define the combustion rig signal analysis 
system the pressure signal measured in the duct, p(t)s due to the combustion 
process, n(t), is modified by rig system operation, g(t). The rig system 
operation function represents the results of the echoes from the combustion 
duct ends and the interaction of the echoes and the combustion source. The 
output of the system can be written as the following convolution integral: 


p{t) 



nfr) g (t 


T 


) dr 


Taking the Fourier transform of Eq. (1) yields 
F |p(t)} = N(f)G(f) 


( 1 ) 


( 2 ) 


Where N(f) is the combustion process spectrum and G(f) is the rig system 
operation transfer function. The pressure spectrum is then given by 


P(f) = 


N(f) 


6(f) 


Taking the logarithm of Eq. (3) yields 


log P(f) = log 


N(f) 


2 + 


log 


6(f)! 


(3) 


(4) 
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The inverse Fourier transform of Eq. (4) yields the power cepstrum 


qUAUTV 


C(q) 



log P(f) 



log 


N(f)f 




G(f) 



( 5 ) 


Provided the combustion source cepstrum and the rig system cepstrum occupy 
different quefrency ranges, the effects of the two systems are separable. 
Consequently, for this condition the individual contributions of each power 
cepstrum can be separated by liftering in the quefrency domain. Thus the 
smoothed cepstrum is 


C 


(q) - 



log 


N(f) 



( 6 ) 


The smoothed spectrum is obtained by Fourier transforming the smoothed 
cepstrum. Thus 


log p’(f) = F|c'(q) 


(7) 


The procedure used to compute the Fourier transform is presented In 
Ref. 6. 


EXPERIMENTAL INVESTIGATION 

The combustion process in a turbofan engine is a source of far-field 
noise. At the NASA Lewis Research Center an extensive program was conducted 
to determine the source and characteristics of combustion noise and its 
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propagation using aircraft gas turbine engines, combustor test facilities 
and ducted combustion rigs (Refs. 8 to 14). Jet A fuel was used in these 
tests. In order to determine if using a gaseous fuel changed the combustion 
noise source or combustion noise propagation, a program was initiated using 
gaseous hydrogen and Jet A for the fuels in a ducted combustion rig. 

The experimental apparatus is shown schematically in Fig. 1. The com- 
bustor rig is shown in Fig. la. The combustor section consists of a 0.2m 
diameter by 0.77m long J-47 burner can placed concentrically in a 0.3m diam- 
eter by 0.77m long flow duct. The combustor section is followed by a 0.203m 
diameter by 6.44m long flow duct. The cut-off frequency of the long duct is 
above 1500 Hz. The Jet A and hydrogen fuel injection schematic is shown in 
Fig. lb. Note, that the hydrogen is injected at the same point in the J-47 
burner can as the Jet A. 

The dynamic pressure measurements were made simultaneously in the duct 
at the four locations shown in Fig. 1 (one just upstream of the combustor 
and three in the long duct at 0.3m, 0.7m, and 6.14m from the long duct in- 
let) and at 3 far-field microphone locations (30, 90, and 120 degrees). 
Typical results from the microphone in the duct at 6.14 m only are presented 
in this paper. 

The internal pressure transducers used were conventional 0.635m micro- 
phones with pressure response cartridges. In order to avoid direct exposure 
of the microphones to the combustion gases, they were mounted outside the 
ducted combustion rig and the fluctuating pressure in the rig was trans- 
mitted to the transducers by means of a "semi-infinite" acoustic wave 
guide. Details on these probes are given in Ref. 8. 

The signals from the four internal probes and three far-field micro- 
phones were FM recorded on magnetic tape in 2 minute record lengths for 
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later processing. The internal probes and far-field microphones were cali- 
brated with a piston phone prior to each day's running. The constant band- 
width spectra given in this paper were obtained by off-line processing of 
the taped data on a two-channel fast Fourier transform digital signal proc- 
essor with a built-in a-d converter and 120dB/octave anti-aliasing filters. 
The processor was capable of direct computation of up to 4096 ensemble aver- 
ages of a 1024 point forward or inverse Fourier transform to yield either 
time-domain or frequency domain information. 

All tests were conducted in an outdoor acoustic arena. The facility is 
shown in Fig. 2. 

The results discussed here were obtained at nominal mass flows of 0.6, 
1.3, and 2.1 kg/s and at nominal temperatures of 760 K, 920 K, and 1090 K 
for hydrogen and Jet. A. The operating conditions are summarized in 
Tables I - III. 


RESULTS 

Measured Spectra 

Typical logarithmic pressure spectra measured in the combustion duct 
near the exit are shown in Figs. 3 and 4 for Jet A and hydrogen fuel, re- 
spectively. The spectrum for Jet A was measured at a combustor exit temper- 
ature of 1093 K at an air mass flowrate of 0.55 kg/s. The spectrum for hy- 
drogen was measured at a combustor exit temperature of 1036 K at an air mass 
flowrate of 0.62 kg/s. Note that the ripples in the hyvVogen spectra have a 
much larger amplitude than those in the Jet A spectrum. Similar variations 
in spectra were observed at the other test conditions covered in the study. 
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In each case the spectrum consists of a sequence of harmonic peaks and dips 


superimposed on an underlying broad-band structure. The spectra have a 
rippled appearance. The peaks of the ripples occur approximately at 
frequencies f^ = (n •♦•l/2)c/2L. Moreover^ the peaks show a regular 
variation of amplitude with frequency. The pattern which is most common is 
as follows: peaks are small at lower frequencies; at intermediate 
frequencies an amplitude maximum is reached; finally, at high frequencies 
the ripple amplitude is small again. 

Cepstra 

Cepstra calculated from the spectra shown in Figs. 3 and 4 are shown in 
Figs. 5 and 6. All the cepstra examined for a range of test conditions have 
a fundamental pattern consisting of a large narrow negative spike near a 
quefrency of 0.03s with two adjacent positive spikes located on either side 
(see Fig. 5). The large excursion near 0.03s corresponds to the reciprocal 
of the spacing between the harmonics of 30 Hz shown in Figs. 3 and 4. Also 
shown in Fig. 5 is a set of smaller excursions near 0.06s which are first 
rahmonics. This pattern is related to the way a plane wave propagates in 
the duct. The cepstra represent the interaction of the following elements: 
amplification in the combustion region, wave reflection at inlet and exit of 
the combustion system, the position of the microphone, the length of the 
duct, and the flow and sound velocities. A mathematical model of the system 
should yield not only autoand cross-spectra but also the cepstra. Such a 
model should jrovide a complete explanation of the fundamental pattern. 

The large narrow negative spike occurs in the fundamental pattern at 
the time (quefrency) it takes a signal to travel up and down the duct 
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dx 
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( 8 ) 


^0 


c(x) 



Furtherinore, adjacent to this negative spike are two positive spikes at 
times (quefrencies) 


= _2_ikil 


c(x) 


X = L 


( 9 ) 



and 


= 2(L-fc) (10) 



where n is the distance between the infinite tube microphone and the 
duct exit. 

The spike at quefrency q^^ is narrow while the spike at quefrency 
(^2 is broad. Also, each cepstrum has an inverse of this fundamental pat- 
tern at the second rahmonic quefrencies Zq^ , 2q^ , and 2q2. More- 
over, in some cases the fundamental pattern is repeated at Sq^ , 3q^ , 
and 3q2» 

The following simple analysis is used to correlate the quefrency of the 
major negative peak in the cepstra with the airflow rate, duct inlet temper- 
ature and duct exit temperature. Because the temperature distribution was 
not measured, it is assumed, for convenience, to be exponential. The duct 
temperature variation with position is assumed to be 
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T(x) . T(x) ^ e ■“(''•^•°^>0<x<7.21m 

where 


( 11 ) 



(12) 


The I'sentroplc speed of sound Is given by 

Cj,(x) = 19.61 |/r(7r (13) 

where the gas specific heat ratio used is 1.4 and is assumed to be indepen- 
dent of X. The flow velocity is assumed to be given by 

v(x) = W ^ T(x) (14) 

A p (MW) 

0 

where the gas molecular weight is also assumed to be independent of x. 

Substituting Eqs. (11) to (14) into Eq. (8), using the measured temper- 
ature and air flows shown in Table III, and numerically integrating yields 
the theoretical echo travel time, (^ q)j i shown in Table III. The 
measured echo travel times corresponding to the quefrency of the major nega- 
tive spike, (qQ)f^ » also shown in Table III. The variance of the 
theoretical echo travel time from the measured echo travel times was about 
20 (ms) . However, the differences between the measured and theoretical 
values showed systematic trends which could be empirically correlated with 
mass flow and temperature using 
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* 1.07m 


(15) 


and the following expcvnents 


e = 0.505 

$ =: -4.53x10-3 
6 = 0.106 


p 

The variance with this equation was reduced to about 1 (ms) ‘. The magni- 
tude of the exponents indicates that a small correction in the temperature 
distribution function might improve the model while changes related to the 
mass flow would make it less accurate. These results show that the cepstrp 
pr'i'^rdure is useful in the study of models for sound propagation in the com- 
bustion duct. 

Furthermore, these results indicate that the sound propagated at the 
normal isentropic speed in the hydrogen and the Jet A fueled tests conducted 
using this facility. However, the echo time delay measurement does not 
provide information on the variation of the sound speed with frequency. An 
analysis of cross spectra, which is beyond the scope of this paper, would 
provide more definitive information on the variation of the speed of sound 
with frequency. 

The fundamental cepstrum pattern (large narrow negative spike with ad- 
jacent positive spikes) appears since part of the measured signal consists 
of three signals which can be considered to have started at the microphone 
and traveled distances 2(L - a) , 2L and 2(L + Ji) before returning 
to the microphone. Each of these signals underwent different reflections 
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along the path it traveled. The inlet of the combustion system acts like a 
closed tube so that the reflection factor at the inlet is R = 1. The exit 
of the duct is assumed to have a R = 1 reflection factor. The spike at q 
is negative because a signal traveling a distance 2L then has a R = 1 
and a R = -1 reflection. Consequently, after the echo delay time the 
direct signal and the echo are 180 degrees out of phase and they will 
combine destructively. However, a signal traveling a distance 2(L - n) 
has only an R = 1 reflection. Thus, the echo undergoes no phase shift. 
While a signal traveling a distance 2(L + a) has an R = 1 reflection 
and two R = -1 reflections and undergoes a 360 degree phase shift. For 
these cases after the corresponding echo delay times the direct signal and a 
copy of the direct signal then will combine constructively. 

Note that the signal that travels a distance 2(L + a) and had an 
R = 1 reflection and two R = -1 reflections is the one responsible for 
the higher quefrency positive broad spike in the fundamental pattern. 
Therefore, for this signal path more dispersion in travel times is present 
than along the signal paths which have one or no reflections from the duct 
exit. 

Most of the hydrogen cepstra show many additional spikes in addition to 
the fundamental pattern (Fig. 6). However, at the higher air flow rates the 
hydrogen cepstra resemble the Jet A cepstra and only the fundamental pattern 
is observed. 


Smoothed Spectra 

Smoothed spectra produced by littering and inverse Fourier transforming 
the cepstra presented in Figs. 5 and 6 are shown in Figs. 7 and 8. For Jet 
A combustion with a combustor exit temperature of 1093 K at a mass flow rate 
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of 0.55 kg/s the duct exit smoothed spectrum is shown in Fig. 7. The 
spectrum has a flat appearance. The smoothed hydrogen spectrum measured at 
a combustor exit temperature of 1036 K at an air mass flow rate of 0.62 kg/s 
is presented in Fig. 8. The hydrogen spectrum has a peaked appearance. 

In order to obtain these smooth spectra, the following liftering pro- 
cedure was used to preserve the broad-band spectral shape while removing 
most of the harmonic ripple in the measured spectrum. Since most of the 
broad bank spectral information is present in the cepstra at quefrencies 
below 0.0046 s, at quefrencies greater than 0.0046 s excursions greater than 
+10 and -10 were set to zero. Then new values were obtained at quefrencies 
which are zero by linear interpolation. 

This liftering procedure is somewhat arbitrary. For example, it is 
possible to remove excursions greater than +5 and -5 or +1 and -1. Somewhat 
smoother spectra are obtained by removing excursions greater than +5 and 
-5. However, removing excursions greater than +1 and -1 may produce smooth 
but distorted spectra if the corresponding cepstrum has many large rahmonics. 

The degree of smoothness and the acceptability of the spectrum after 
liftering are subjective decisions. However, for a wide range of choices 
the resulting overall broad-band power levels after lifterinng are in good 
agreement. For a typical spectrum having tones where the maximum level 
change between a cancellation and a reinforcement is 10 dB, the resulting 
overall power level before smoothing was 142.02 dB, Liftering by removing 
excursions of +10 and -10, +5 and -5, and +1 and -1 in the cepstrum yields 
overall broad-band power levels of 140.72, 140.60, and 140.57, respec- 
tively. As another example, for a typical spectrum having tones where the 
maximum level change between a cancellation and a reinforcement is 30 dB, 
the resulting overall power level was 146.92 dB. Liftering by removing 
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excursions of +10 and <-10, +5 and -5, and +1 and -1 in the cepstrum yields 
overall broad-band power levels of 133.48, 133.53, and 134 db, re- 
spectively. Consequently, for the cases studied the scatter in the overall 
broad-band power levels obtained from smoothed spectra by using any of these 
sets of limiting values in the liftering procedure is less than 1 dB. 

The significant differences between the Jet A and hydrogen spectra are 
not due to propagation since (1) the fundamental cepstrum pattern (large 
narrow negative spike with adjacent positive spikes) occurs in both sets of 
data and (2) the sound propagates at the normal isentropic speed in both Jet 
A and hydrogen. Consequently, in order to understand why the measured vari- 
ation of the spectra with operating condition occurs, the combustion process 
and the combustion noise must be studied simultaneously. 

Correlation with Heat Release Rate 

The overall acoustic broad-band power levels (OAPWL) were obtained from 

the smoothed duct spectra. In calculating the overall pressure levels the 

fluctuating pressures are assumed to be acoustic. The sound power level 
-13 

(re. 10 watts) is calculated assuming plane wave propagation of the 
pressure perturbation. Various methods for relating the total pressure 
level to the steady state operating conditions have been described in 
Refs. 15 and 16. 

The correlation of overall acoustic broad-band power level and heat 
release rate squared is shown in Fig. 19 for Jet A and hydrogen where 
Q = W c aT. It is apparent that significant scatter exists in the 

H 

acoustic power level data. However, for this system operating at 
atmospheric pressure and low heat release rates, the overall power level 
dependence on the square of the heat release rate is consistent with previ- 
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ous data. For example, this same dependence was also found In a study of 
core noise turbofan data from can combustors (Ref. 16) and annular combus- 
tors (Ref. 9) for the low heat release rates. Furthermore, the hydrogen 
overall acoustic broad-band power level for low heat release rates is about 
3 dB greater than the Jet A overall acoustic broad-band power level at the 
same heat release rate. 


CONCLUSIONS 


Based on the application of cepstral techniques to the study of the 

pressure auto spectra measured near the exit of a long duct connected to a 

J-a.7 combustor using both hydrogen and Jet A fuels it is concluded that: 

(1) The cepstra have a fundamental pattern (large narrow negative spike 
with adjacent positive spikes) due to signal propagation along three 
different paths. 

(2) The fundamental pattern in the cepstra can be used to extract informa- 
tion on sound propagation in the duct. 

(3) The quefrency of the large negative spike in the measured cepstra is 
consistent with the isentropic propagation of sound in both the Hydro- 
gen and Jet A. 

(4) Organ pipe tones can be removed from measured spectra by a liftering 
technique. 

(5) The scatter in the overall broad-band power levels obtained from spec- 
tra smoothed by using three liftering techniques is less than 1 dB. 

(6) Differences between the Jet A and hydrogen spectra for the same test 
condition are due to the combustion source region and not propagation 
effects. 
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(7) The variation with operating condition of the overall acoustic broad- 
band power level for tests using both hydrogen and Jet A is consistent 
with previous data that found it was proportional to the square of the 
heat release. 

(8) The hydrogen overall acoustic broad-band power level for low heat re- 
lease rates is greater than the Jet A overall acoustic broad-band power 
level at the same heat release rate. 
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TABLE I. - JET A TEST CONDITIONS 


Test point 

w, 

kg/s 

COMBUSTOR 

Q, 

MW 

P, 

N/m^ 

OASPL, 

dB 

^AIR, 

'K 

^EXIT, 

•k 

61 

0.58 

278 

761 

0.31 

1.01 

124.80 

59 

.57 

279 

916 

.42 

1.02 

125.50 

37 

.55 

278 

1093 

.54 

1.01 

126.40 

44 

1.34 

278 

759 

.66 

1.08 

132.90 

46 

1.27 

278 

924 

.95 

1.09 

134.40 

48 

1.2G 

279 

1090 

1.23 

1.10 

136.70 

54 

2.11 

278 

756 

1.12 

1.25 

141.60 

56 

2.12 

279 

925 

1.58 

1.27 

143.10 

58 

2.13 

279 

1092 

2.07 

1.28 

143.10 
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TABLE II. - HYDROGEN TEST CONDITIONS 


Test point 

w, 

kg/s 

COMBUSTOR 

Q, 

MW 

P, 

N/M^xlO^ 

PWL, 

dB 

’'’air, 

‘k 

"^EXIT, 

“K 

30 

0.57 

276 

755 

0.30 

1.01 

128.70 

8 

.63 

276 

923 

.48 

,990 

133.50 

9 

.62 

276 

1036 

.57 

.991 

136.40 

3 

1.29 

276 

755 

.69 

1.00 

139.50 

1 

1.26 

276 

924 

.96 

1.00 

140.03 

6 

1.30 

277 

1092 

1.28 

1.01 

138.30 

22 

2.14 

276 

756 

1.15 

1.04 

138.60 

24 

2.15 

277 

922 

1.62 

1.04 

141.10 

21 

2.13 

276 

1098 

2.12 

1.06 



142.90 


















TABLE III. - RESULTS OF ECHO TRAVEL TIME ANALYSIS. 
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Figure 2. - NASA Lewis combustion acoustics facility. 



Figure 3. - Duct exit spectra (bandwidth • 2 
Hz) for Jet A combustion at a combustor 
exit temperature of 1093 K at a mass flow 
of 0.55 kg/s. 
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Figure 4 - Duct exit spectra (bandwidth ■ 
2 Hz) for hydrogen combustion with a 
combustor exit temperature of 1036 K 
at a mass flow rate of 0.62 kg/s. 
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Figure 5. - Duct exit cepstra for Jet A combusti^jn 
with a combustor exit temperature of 1093 K at 
a mass flow rate of 0. 55 kg/s. 
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Figure 6. - Duct exit cepstra for hydrogen combus- 
tion with a combustor exit temper? tu re of 1036 K 
at a mass flow rate of 0.62 kg/s. 
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Figure 7. - Smoothed duct exit spectrum 
for Jet A combustion with a combustor 
exit temperature of 1093 K at a mass 
flow rate of 0. 55 kg/s. 
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Figure 8. - Smoothed duct exit spectra for 
hydrogen combustion with a combustor 
exit temperature of 1036 K at a mass 
flow rate of 0.62 kg/s. 
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Figure 9. ~ Correlation of combustion noise and square of heat release 
rate, Q (Q in megawats). 



